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0\ . ABSTRACT 

rG ' We present here initial results on the X-ray Cluster Luminosity Function (XCLF) from the 

Bright SHARC sample of distant X-ray clusters of galaxies. This sample is 97% complete in its 



optical identifications and contains 12 X-ray luminous clusters in the redshift range 0.3 < z < 0.83 
(median z = 0.42) and 1.1 < L x < 8.5 x 10 44 ergs _1 [0.5-2.0 keV]. We present a preliminary 
selection function for the Bright SHARC based on Monte Carlo simulations. Using this selection 
function, we have computed the Bright SHARC XCLF and find it to be fully consistent with a 
non-evolving XCLF to L x ~ 5 x 10 44 ergs _1 and z ~ 0.7. 

At L x > 5 x 10 44 ergs -1 , we find evidence for a deficit of clusters compared to that expected 
from a non-evolving XCLF. We detect only one such cluster in the redshift range 0.3 < z < 0.7 
when we would expect 4.9 clusters based on the local XCLF of De Grandi et al. (1999). The 
statistical significance of this deficit is ~ 96%. To increase the statistical significance of this 
possible deficit, we have combined the Bright SHARC and the 160 deg 2 survey of Vikhlinin et 
al (1998a). This joint survey covers ~ 260 deg 2 and contains only one confirmed 0.3 < z < 0.7, 
L x > 5 x 10 44 ergs _1 cluster, while we would expect 7.6 such clusters based on the local XCLF 
(De Grandi et al. 1999). The statistical significance of the deficit in this joint survey increases 
to 99.5%. These results remain preliminary because of incompletenesses in the optical follow-up 
and uncertainties in the local XCLF. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: evolution - 
surveys — X-rays: general 



1 Based on data obtained at the Kitt Peak National Observatory, European Southern Observatory, the Canada-France-Hawaii 
Telescope and Apache Point Observatory 



1. Introduction 

The observed evolution of the space density of clusters of galaxies provides a powerful constraint on 
the underlying cosmological model. Many authors have demonstrated - both analytically and numerically - 
that the expected abundance of clusters, as a function of cosmic epoch, is a sensitive test of the mean mass 
density of the universe (fi m ) and the type of dark matter (Press & Schechter 1974; Lacey & Cole 1993, 1994; 
Oukbir & Blanchard 1992, 1997; Henry 1997; Eke et al. 1996, 1998; Viana & Liddle 1996, 1999; Bryan & 
Norman 1998; Reichart et al. 1999; Borgani et al. 1999). 

Measurements of the evolution of the cluster abundance have made significant progress over the past 
decade. For example, in their seminal work, Gioia et al. (1990) and Henry et al. (1992) computed the 
luminosity function of X-ray clusters extracted from the Einstein Extended Medium Sensitivity Survey 
(EMSS) and concluded that the X-ray Cluster Luminosity Function (XCLF) evolved rapidly over the redshift 
range of 0.14 < z < 0.6. 

The launch of the ROSAT satellite heralded a new era of investigation into the XCLF. The ROSAT 
All-Sky Survey (RASS) has provided new determinations of the local XCLF and has demonstrated that 
there is little observed evolution in the XCLF out toz~ 0.3 (Ebeling et al. 1997; De Grandi et al. 1999) in 
agreement with the earlier work of Kowalski et al. (1984). In addition, the ROSAT satellite has supported 
several investigations of the distant X-ray cluster population (RIXOS, Castander et al. 1995; SHARC, Burke 
et al. 1997, Romer et al. 1999; RDCS, Rosati et al. 1998; WARPS, Jones et al. 1998; Vikhlinin et al. 1998a; 
NEP, Henry et al. 1998). Initially, such investigations reported a deficit of high redshift, low luminosity 
clusters consistent with the original EMSS result (Castander et al. 1995). However, over the last few years, 
there has been a growing consensus for a non-evolving XCLF. First, Nichol ct al. (1997) re- examined the 
EMSS cluster sample and determined that the statistical evidence for evolution of the EMSS XCLF had 
decreased in light of new ROSAT data. Second, several authors have now conclusively shown that the XCLF 
does not evolve out toz~ 0.7 for cluster luminosities of L x < 3 x 10 44 ergs _1 (Collins et al. 1997; Burke et 
al. 1997; Rosati et al. 1998; Jones et al. 1998). 

Above L x = 3x 10 44 ergs _1 , recent work has indicated that the XCLF may evolve rapidly in agreement 
with the original claim of Gioia et al. (1990). Reichart et al. (1999) highlighted a deficit of luminous 
(L x > 5 x 10 44 ergs _1 ) EMSS clusters at z > 0.4 i.e. the EMSS survey has both the sensitivity and area 
to find such clusters but does not detect them. Moreover, Vikhlinin et al. (1998b) has recently reported 
evidence for a deficit of luminous clusters at z > 0.3 based on the 160 deg 2 ROSAT survey (Vikhlinin et al. 
1998a). 

In this paper, we report on the first determination of the bright end of the XCLF that is independent of 
the EMSS. In sections || & y, we outline the Bright SHARC sample of clusters used herein and its selection 
function. In sections [| & [^, we present the derivation of the XCLF and discuss its implications. Throughout 
this paper, we use H = 50km s _1 Mpc and q = ^ to be consistent with other work in this field. All quoted 
luminosities are in the hard ROSAT passband [0.5 — » 2.0 keV] and are aperture and k-corrected (see Romer 
et al. 1999 for details). 



2. The Bright SHARC Sample 

The details of the construction of the Bright SHARC survey are presented in Romer et al. (1999). The 
Bright SHARC was constructed from 460 deep (T cxp > 10 ksecs), high galactic latitude (|6| > 20°), ROSAT 



PSPC pointings which cover a unique area of 178.6 deg 2 . Using a wavelet-based detection algorithm, 10277 
X-ray sources were detected in these pointings of which 374 were measured to be significantly extended 
(> 3er; see Nichol et al. 1997) relative to the ROSAT PSPC point-spread function. The Bright SHARC 
represents the brightest 94 of these 374 extended cluster candidates above a ROSAT count rate of 0.0116 
cntss -1 . This corresponds to a flux limit of ~ 1.4 x 10~ 13 ergs -1 cm 2 [0.5-2.0 keV] for the average neutral 
hydrogen column density of the Bright SHARC and a cluster temperature of 6 keV. 

Over the past two years, we have optically identified the most likely X-ray emitter for 91 of these 94 
Bright SHARC cluster candidates and have discovered 37 clusters, 3 groups of galaxies and 9 nearby galaxies 
(the remainder are blends of X-ray sources e.g. AGNs & stars; see Romer et al. 1999). We find 12 clusters 
in the range 0.3 < z < 0.83 (median redshift of z = 0.42) and have independently detected cluster RXJ0152- 
7363 (z = 0.83 based on 3 galaxy redshifts obtained at the CFHT) which is one of the most luminous, high 
redshift X-ray clusters ever detected (see Romer et al. 1999). This cluster has also been detected by the 
WARPS and RDCS surveys (see Ebeling et al. 1999; Rosati, private communication). 



3. Selection Function 

An important part of any survey is a solid understanding of the selection function i.e. the efficiency 
of finding objects as a function of both cosmological and operational parameters. In the case of the EMSS 
cluster sample, the selection function is somewhat straightforward since the EMSS optically identified all 
sources regardless of their observed X-ray extent. This is not the case for the Bright SHARC and therefore, 
the most direct way of modelling the selection function is through Monte Carlo simulations. The details of 
such simulations are given in Adami et al. (1999) but we present here some initial results. 

The Bright SHARC selection function is obtained by adding artificial clusters to PSPC pointings and 
determining if these clusters would have satisfied the Bright SHARC selection criteria. In this way, we 
can accurately model the effects of the signal-to-noise cut, the extent criteria and blending of sources. To 
date, we have only simulated artificial clusters at four different intrinsic luminosities, L x = 1, 2, 5 & 10 
(x 10 44 ergs _1 i.e. equally spaced in log-space), and covering a redshift range of 0.3 < z < 1 in steps of 
S z = 0.05. These artificial clusters were constructed assuming an isothermal King profile with r c — 250 kpc 
and (3 = |. The effects of changing the cluster profile and its parameters are explored in detail in Adami et 
al. (1999). 

For each combination of L x and z, we added, one at a time, 10 artificial clusters to each of 10 PSPC 
pointings (randomly chosen from all pointings available in the SHARC survey) . The positions of the clusters 
in these pointings were chosen at random but recorded for later use when computing the area surveyed 
(see below). For each artificial cluster, we computed the expected number of ROSAT PSPC counts for its 
given luminosity, redshift and the exposure time of the pointing. We then took a Poisson deviate about the 
expected number of counts and distributed the counts at random assuming a redshifted King profile. The 
pointing was then processed in the same fashion as the real data (see Romer et al. 1999; Nichol et al. 1997) 
thus allowing us to determine if the cluster would have been selected for the Bright SHARC. 

In Figure El we show the results of these initial simulations. The effective area of the Bright SHARC was 
computed by splitting the PSPC field-of-view into 4 annulii (2/5 -» 6/25, 6/25 -> 11/25, 11/5 -> 16/25 
and 16/5 — > 22/5; the central part of the PSPC was excluded) and multiplying the area in each of these 
annulii by the measured success rate of detecting clusters in these same annulii (as a function of the input 
cluster redshift and luminosity) . We then summed the effective area in these annulii over all pointings used - 
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Fig. 1. — The effective area of the Bright SHARC sample of clusters as a function of cluster luminosity and 
redshift. These curves were constructed from the initial results of our Monte Carlo simulations (see Adami 
et al. 1999). The four curves represent the four different input luminosities (in units of 10 44 ergs _1 ). The 
simulations were performed over a redshift range of 0.3 < z < 1 in steps of S z = 0.05. 

again as a function of luminosity and redshift - to provide the total effective area sampled by the simulations. 
Finally, we re-scaled the results to obtain the expected area for all 460 PSPC pointings used in the Bright 
SHARC. 



4. Luminosity Function 

The luminosity function of the Bright SHARC was determined using the 1/V a methodology outlined 
in Avni & Bahcall (1980), Henry et al. (1992) and Nichol et al. (1997), where V a is the available sample 
volume for any given cluster in the survey. Using the selection function presented in Figure |lj, V a can be 
computed for a cluster of luminosity L x using 



V„ 



Zhigh 



fl(L K , z) V{z) dz, 



(1) 



where zi ow and Zhigh are the lower and upper bounds of the redshift shell of interest, V(z) is the volume per 
unit solid angle for that redshift shell and 0(L X , z) is the effective area of the Bright SHARC from Figure [I]. 
In practice, the integral in Eqn. [fl is replaced by a sum over the discrete values of ft(L x , z) obtained from 
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Fig. 2 — The XCLF for the Bright SHARC survey (solid stars) and the Southern SHARC (solid circles; 
taken from Table 1 of Burke ct al. 1997). The y-axis was computed using Equation g with AL values quoted 
in Table [l] and has units of Mpc~ 3 (10 44 ergs -1 ) -1 . These data points span a redshift range of 0.3 < z < 0.7. 
The hashed region represents the local XCLF as measured by De Grandi et al. (1999) and Ebeling et al. 
(1997); we have included the ±1ct errors as given by these authors. The unfilled stars are taken from Nichol 
et al. (1997) and are the EMSS luminosity function at 0.3 < z < 0.6 including new data from the ROSAT 
pointing archive. The solid triangle is the XCLF expected for RXJ0152-7363 and its value is presented in 
Table [jL All the error bars are computed using the Poisson errors of Gehrels (1986) 

the simulations. Linear interpolation was used where necessary to obtain finer resolution in both luminosity 
and redshift space. 

The luminosity function was derived by summing the V a values for all clusters in the Bright SHARC as 
a function of luminosity i. e. 
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where AL is the width of the luminosity bins and N is the number of clusters in that luminosity bin (Table |lj). 
For the results presented here, we have restricted ourselves to the luminosity range 10 44 < L x < 10 45 ergs _1 . 

In Figure |, we present the Bright SHARC XCLF and compare it to the Southern SHARC XCLF (Burke 
ct al. 1997) and measurements of the local XCLF (Ebeling et al. 1997; De Grandi et al. 1999). We provide, 



Table 1: The measured space density of the Bright SHARC sample of clusters 
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44.23 




-6.74 


1.4 




3 
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in Table [jj the data points displayed in Figure |] together with the redshift and luminosity ranges studied. 
We also provide the number of clusters in each bin. We have not performed a parametric fit to the data 
because of the limited dynamic range in luminosity available from our present simulations. 



5. Discussion 

Figure |] demonstrates that the high redshift XCLF does not evolve below L x = 5 x 10 44 erg s _1 , or below 
L* in the XCLF (L* = 5.7±J;^ x lO^ci-gs^ 1 from Ebeling et al. 1997). Using a Kolmogorov-Smirnov (KS) 
test similar to that discussed by De Grandi et al. (1999), we find that the Bright SHARC unbinned data 
is fully consistent with the low redshift XCLF over the luminosity range 1 < L x < 5 x 10 44 crgs _1 (we find 
a KS probability of 0.32). This result is consistent with previous work (Nichol et al. 1997; Burke et al. 
1997; Rosati et al. 1998; Jones et al. 1998) but pushes the evidence for a non-evolving XCLF to the highest 
luminosities presently reached by ROSAT data. This is illustrated by the fact that we find 6 distant clusters 
in the luminosity range 3 < L x < 5 x 10 44 ergs _1 (Table fll) which is more than any other ROSAT archival 
survey. The issue therefore becomes the degree of observed evolution in the XCLF above L*. 

We have investigated evolution in the XCLF above L* in two separate ways. First, we have one very 
luminous high redshift cluster in the Bright SHARC - RXJ0152.7-1357 (z = 0.83; L x = 8.5 x lO^ergs- 1 ) 

- which probes the XCLF above L*. The implied position of this cluster in the XCLF is given in Figure @ 
and Table p] (we used a redshift range of 0.3 < z < 1 and luminosity range of 7 x 10 44 < L x < 10 45 ergs" 1 
when computing the volume sampled by this cluster) . As can be seen in Figure @, the implied space density 
of RXJ0152-1357 agrees with the local XCLF and may be evidence for a non-evolving XCLF above L*. 
However, we must remain cautious since RXJ0152. 7-1357 has a complex X-ray morphology indicative of 
an on-going merger which may have enhanced its luminosity (see Ebeling et al. 1999). Such disturbed or 
non-spherical morphologies (both in the X-rays and optical) appear to be common at these high redshifts 

- e.g. MS1054.4-0321 at z = 0.823 (Donahue et al. 1999) and RXJ1716. 6+6708 at z = 0.813 (Henry et al. 
1998; Gioia et al. 1999) - and may indicate that we are witnessing the epoch of massive cluster formation. 

The second path of investigation is to repeat the analysis of Collins et al. (1997) and Vikhlinin et 
al. (1998b) and compute the number of expected Bright SHARC clusters at these bright luminosities 
assuming a non-evolving XCLF. In the luminosity range 5 x 10 44 < L x < 10 45 ergs _1 and redshift range 
0.3 < z < 0.7, we would predict 4.9 clusters based on the De Grandi et al. (1999) XCLF (or 3.5 clusters 
using the Ebeling et al 1997 XCLF). At present, the Bright SHARC contains only one confirmed cluster in 
this range, RXJ1120. 1+4318 at z = 0.60 and L x = 5.03 x 10 44 ergs _1 . The Poisson statistical significance 
(Gehrels 1986) of this observed deficit is ~ 96% (or ~ 90% for Ebeling et al. 1997). 



One way to increase the statistical significance of any possible deficit of high redshift X-ray luminous 
clusters is to combine the Bright SHARC and 160deg 2 survey of Vikhlinin et al. (1998a) as both surveys 
should have similar selection functions. We have determined that 201 ROSAT PSPC pointings are in common 
between the two surveys (see Romer et al. 1999), or ~ 44% of the area of the Bright SHARC. Vikhlinin et al. 
(1998b) have noted that their survey probably contains only 2 candidate clusters above L x = 3x 10 44 ergs _1 
and z > 0.3 (based on photometric redshift estimates), while they would have expected 9 clusters above this 
redshift and luminosity limit. In fact, Bright SHARC spectroscopy of one of these two candidate clusters - 
RXJ1641+8232 - shows that it is at z — 0.195 (based on 3 galaxy redshifts). The other candidate cluster - 
RXJ1641+4001 - is not in the Bright SHARC. 

Therefore, the Bright SHARC plus the 160deg 2 survey cover a total area of ~ 260 deg 2 and scaling the 
aforementioned numbers appropriately, we would expect 7.6 (using the De Grandi et al. 1999 XCLF) or 5.5 
(using the Ebeling et al. 1997 XCLF) clusters in the luminosity range 5 x 10 44 < L x < lO^ergs" 1 and 
redshift range 0.3 < z < 0.7. By comparison, this joint survey contains only 1 confirmed 2 ] X-ray luminous, 
high redshift cluster. From Gehrels (1986), the statistical significance of this deficit is now 99.5% or 97.5% 
respectively. We note however that there are still candidates in both surveys that require further optical 
follow-up i.e. one of the three unidentified Bright SHARC candidates mentioned in section could be a 
high redshift cluster. 

The above analyses suffer from small number statistics and incomplete optical follow-up. Moveover, the 
local XCLFs are also affected by small number statistics at these high X-ray luminosities. Therefore, the 
issue of evolution above L* in the XCLF remains unclear. However, we note that four independent surveys of 
distant clusters of galaxies have now seen a potential deficit of X-ray luminous high redshift clusters (SHARC, 
WARPS, EMSS, 160 deg 2 survey). It is worth stressing that some level of XCLF evolution is expected above 
L* at high redshift and the degree of such evolution is a strong indicator of Vt m (see Blanchard & Oukdir 
1992, 1997; Reichart et al. 1999). 

The way to resolve these problems is to construct larger samples of clusters over bigger areas of the 
sky. This is certainly possible for the local XCLF e.g. using the on-going REFLEX survey of clusters being 
constructed from the RASS (see Bohringer et al. 1998). However, for the distant XCLF, it is unlikely that 
a significant amount of further area (>> 260 deg 2 ) can be added to the present surveys using the existing 
ROSAT pointing archive; one may have to wait for the XMM satellites (see Romer 1998). The next major 
improvements to the results presented in this paper will be to obtain a more detailed view of the Bright 
SHARC selection function as well as to complete the optical follow-up of the remaining candidates. Another 
improvement would be the possible combination of all the ROSAT distant surveys (SHARC, 160 deg 2 , RDCS, 
NEP & WARPS) thus maximising the amount of volume sampled at high redshift. 
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